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Abstract. The present understanding of the electronic and magnetic properties of a'-NaV^Os is based 
on the hypothesis of strong charge disproportionation into V 4+ and V 5+ , which is assumed to lead to a 
spin-1/2 Heisenberg chain system. A recent structure analysis shows, however, that the V-ions are in a 
mixed valence state and indistiguishable. We propose an explanation for the insulating state, which is not 
based on charge modulation, and show that strong correlations together with the Heitler-London character 
of the relevant intermediate states naturally lead to antiferromagnetic Heisenberg chains. The interchain 
coupling is weak and frustrated, and its effect on the uniform susceptibility is small. 
Dedicated to J. Zittartz on the occasion of his 60th birthday 

PACS. 75.50.Ee - 75.30.Et - 75.40.Cx - 75.40.Mg 

One-dimensional spin-1/2 Heisenberg antiferromagnets determined by neutron scattering ||, Raman|| and other 
are expected to undergo a structural phase transition into experiments. The transition into the low-temperature dimer- 



| a dimerized phase at low temperature accompanied by the ized structure was confirmed by X-ray scattering^ , N MR|tJ , 

_ i 

opening of a spin gap[Q. This spin-Peierls transition was thermal-expansion || yet the detailed deformation pattern 

first observed in organic systems but has found con- is still unknown. 

siderable experimental attention after its recent discovery Based on an early structure determination for a'-NaV 2 5 

in CuGe0 3 (T sp = 14K) j| . The a'-phase of NaV 2 5 ap- |j| i the current picture for the origin of the one-dimensional 

pears to be the second inorganic compound where a simi- magnetic properties rests on the assumption of charge dis- 

lar transition was observed with an even higher transition commensuration into V 4+ and V 5+ chains [§,|,0,|lO|,0. In 

temperature T sp = 3AK 0. The size of the spin gap was Fig.l the VI and V2 chains in b-direction would corre- 
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spond to these different charged vanadium chains, fn this 
picture the V 4+ chain would correspond to a spin- 1/2 
Heisenberg chain. A problem with this picture is, however, 
the missing physical argument for such a strong charge 
modulation. 

In fact a recent new structural analysis of a'-NaA^Os 
by von Schnering et a/ fl2]| clearly shows that all vanadiums 
are equivalent and in a mixed valent state. Furthermore 
the charge-ordered state at room temperature was found 
to be in conflict with Raman scattering experiments ||l3f . 
Since there is one d-electron per two V ions an explanation 
for the insulating properties must be found which does not 
rest on charge discommensuration. 

The aim of this work is to show how in such an in- 
termediate valence situation correlation effects can lead 
to a quasi-lD Heisenberg antiferromagnet. An important 
structural information [fl2"f , which is crucial for our analy- 
sis, are the V-0 distances d in the a-b plane (Fig.l): d(Vl- 
03)=d(V2-03)=1.825A, d(Vl-02||b)=1.916l and d(Vl- 
02||a)=1.986l. This implies that the V1-03-V2 bond is 
much shorter than the V1-02-V1 bond in b-direction. This 
additional bonding effect can be attributed to the d xy - 
electrons. From Fig. 1 we see that the d^-orbitals have a 
direct overlap in one direction, e.g. V2-V1', which is how- 
ever quite small t xy = 0.3 eV|H. Yet a d^-electron can 
also hop between VI and V2 via a double exchange like 
process. In such a process first an electron hops from the 
occupied 03 Py-orbital to V2 with an excitation energy 
At y = e(d xy ) — e(p y ), and in a second step the Vl-electron 
annihilates the oxygen hole. The matrixelement for this 
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Fig. 1. Orbital structure of a'-NaV^Os in the a-b plane. The 
lowest d-orbitals with d^-symmetry are occupied by one elec- 
tron per two V-atoms. The O p-orbitals at the corners of square 
pyramids (where only the basal plane is indicated) are occu- 
pied, with energies ranging from -3 eV down to -7 eV. Solid 
(dashed) squares indicate the downward (upward) orientation 
of the pyramids (for structural details see Refs.[12] and [4]). 
The d xy orbitals have a direct overlap t xy along the V2-V1'- 
V2-V1' zig-zag chain. The largest hopping matrix element t a 
is, however, via double exchange interaction involving an 03 
p y orbital. The resulting bonding is manifested in the struc- 
ture by the shortest planar V-0 distance along the V1-03-V2 
bond. 

process is t a = t 2 d /Ae y . The matrixelement t p d depends 
on the vanadium-oxygen distance d and can be estimated 
with the help of Harrison's Jl5| solid state table as t p d = 
Vpdirtfrl^/md 7 / 2 , where r/ pd7r = 1.36, h 2 /m = 7.62 eVl 2 , 
and Td — 0.98A for vanadium. This gives t p d = 1.2 eV. 
The p-d excitation energies are determined from a LMTO- 
bandstructure calculationfl^j, which yields Ae y = 4 eV 
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and Ae x = 6.5 eV, respectively. Hence t a = 0.35 eV, while 
the corresponding matrixelement t b = 0.15 eV is consid- 
erably smaller because of the larger oxygen- vanadium dis- 
tance and the larger Ae x . We note that the larger value 
for t a is consistent with the stronger bonding of the VI- 
03- V2 bond. 

The Hamiltonian for the d^-electrons may be written 
in terms of creation and density operators d\ and n\ aa = 
d iaa di ar7 , respectively, in the form 



H d = - ^2t a {d\ la d i2 a + H.c.) + t/ d y^ w iaT ni, 



i'-> 1 



(ij)a/3er 



(1) 



where we have introduced a cell structure. Here a cell con- 
tains two vanadium atoms VI and V2, i.e. a = 1 and 
2, and is labeled by a cell index i = (i a ,ib)- The cell- 
Hamiltonian consists of a kinetic energy term t a and the 
local interaction Ua — 4 eV[|l6|]. In Hd we droped a shift 
of the single particle levels which is only relevant for to- 
tal energy considerations Jt7| . The hopping between cells 
is defined by the last term, where 



t b for a = (3, j a = i a , jb = h ± 1 



(2) 



t xy for (3, j a =i a ±l, j b = i b ± \ 
Since the ratios between Ud and the various hopping ma- 
trix elements are quite large, we are confronted with a 
strong correlation problem. 

In the following we shall use a cell-perturbation method 
||l8f based on the above cell structure, where each VTV2 
cell contains in the average a single d-electron. The ad- 
vantage of the cell decomposition is the capability to treat 
the local correlations exactly. The complete set of states 
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Fig. 2. Geometric structure of the effective spin model. A site 
corresponds to a V1-03-V2 cell with a single d-electron. Dark 
symbols indicate the part of the structure shown in Fig. 1. 
The largest antiferromagnetic interaction is Jb, while the frus- 
trated interaction J a between neighboring b-chains is small. 
The one-dimensional magnetic structure is further enhanced 
by the topology of the lattice. The antiferromagnetic interac- 
tion J a of spin 3 with spins 1 and 2 on the neighbor chain is 
frustrated. 

is then labeled by the quantum numbers of the cells, i.e. 
including the number of electrons within a cell. In the 
one electron sector the cell states are simply bonding and 
antibonding states at energies ±t a and corresponding op- 
erators bl a (a| CT ) = ^(4i<t ± 4l<t)> respectively. The low 
energy configurations b\ a b^ al of two electrons in neighbor 
cells i and j are coupled in second order due to the hopping 
tiy The intermediate states have two electrons in one cell. 
The low-energy singlet and triplet states have excitation 



energies AE S = 2t a -±{^fUj+m? a -U d ) and AE t = 2t a , 
respectively. Other singlet states are at much higher en- 
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ergy E > Ud + 2t a . The Heitler-London singlet-triplet 
splitting of the low-energy intermediate states turns out 
to be crucial for the anisotropic nature of the magnetic 
properties. 

The coupling of the b\ a b^ a , configurations may be ex- 
pressed in compact form by the spin- 1/2 Hamiltonian 

H = E[ J y( s ' s J - i^i) - 4( s i s j + !^i)]> ( 3 ) 

<ij) 

1 4. 

where Si = jb[ a T aa >b\, a i defines the spin of bonding elec- 
trons in terms of the vector of Pauli spin matrices r and 
n\ = J2 a b[ a b\,a their density. Here Jjj and Jy denote the 
coupling of neighboring cells in the singlet- and triplet 
channel, respectively. The relative size of the exchange in- 
tegrals is strongly influenced by the geometrical structure 
(Fig. 1 and 2). The coupling of cells in a-direction, i.e. via 
t xy , yields J s a = 2{t xy /2f/AE s and J* = 2{t xy /2f / AE t . 
They differ only because of the singlet-triplet splitting 
E s - t = \(^Uj + 16*| - U d ). The coupling of cells in b- 
direction is = 2t\/AE s and J\ — 0, where the latter 
exchange integral vanishes due to symmetry. This has the 
effect that the total exchange constants in H = ^ ^ijSiSj 
almost cancel along a-direction, i.e. J a — J* — J*, whereas 
along the 'b-chains' Jb — J£ there is no such reduction. 
These estimates are valid if the ratios ^ and with 
*' = {*;,, t xy /2}, are small compared to 1. Numerical es- 
timates for these exchange integrals based on degenerate 
perturbation theory are shown in Fig. 3 as function of the 
vanadium-oxygen (Vl-03) hybridization t p d- For the es- 
timated value t p d — 1.2 eV one finds Jb — 75 meV and 
J a — 13 meV. 
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Fig. 3. AF-exchange coupling between d-electrons in near- 
est neighbor V1-03-V2 bonds as function of t p d\ (a) Jb along 
b-direction (solid line), (b) the small coupling J a = J% — J a 
between different b-chains mediated by t xy ~ 0.3 eV results 
from the cancellation of the contributions and J a from 
the Heitler-London split singlet and triplet intermediate states 
(dashed lines). 

Experimental estimates for the exchange constants are 
usually obtained from the position of the maximum of 
the uniform susceptibility x(T). In the following we study 
the effect of the interchain coupling J a and the thereby 
introduced frustration of the spin-1/2 model (3) using a 
finite temperature diagonalization technique [L^. The ge- 
ometrical structure of the lattice (Fig. 2) is similar to 
the resonating valence bond systems studied by Anderson 
and Fazekas [Q. Results for a two-leg ladder ('railroad- 
trestle') with periodic boundary conditions along a- and b- 
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Fig. 4. Uniform susceptibility x(T) for a two-leg ladder (2x12) Fig. 5. Uniform susceptibility for a 2 x 12 two-leg ladder 

with the structure given in Fig. 2 for different interchain cou- for different interchain coupling strength J a /Jb = 0,0.2,0.4 

pling strength J a /Jb = 0,0.2 and 0.4 and periodic boundary and 1.0 on a large temperature scale. The rapid drop below 

conditions. T/J b < 0.2 is a finite size effect. 

direction and different interchain coupling strength J a /Jb cent study of this effect by Sandvik et aZ.||| shows that 
are shown in Fig. 4 and 5. Since the interchain coupling this may lead to a reduction of T™ ax by 15% for a one- 
is frustrated the change of \{T) is relatively small. The dimensional Heisenberg chain. Combining these two ef- 
maximum of x(T) is at the temperature T™ ax = a x J b , fects one arrives at the estimate a^ ot ~ 0.50. From the 
where a x = 0.8 for a 2 x 12 system in the absence of the susceptibility measurements P,|lO[ for a'-NaV 2 5 one finds 
interchain coupling, i.e. J a /J b = 0. The exact result for T™ ax ~ 350 K, which leads to the estimate J^ xp ~ 700 K. 
the thermodynamic limit recently obtained by Eggert et This is in reasonable agreement with the theoretical value 



aL|2l| using the thermal Bethe ansatz is a x = 0.6. Inter- J b ~ 75 meV derived in this work. In view of the simplicity 



chain coupling leads to a small shift of the maximum by of the model and the rough estimates of the parameters 

about 2% to lower temperatures for the value J a / J b ~ 0.2 this is quite satisfactory. 

estimated above (Fig. 4). We briefly comment on the band picture, i.e. ignoring 

We note that lattice fluctuations are expected to lead the effect of electron correlations. Since double-exchange 

to a further weakening of antiferromagnetism and an ad- is not contained in this scheme, the only large hopping 

ditional shift of the maximum of x(^) to lower T. A re- matrix element is t xy . This leads to two degenerate one- 
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dimensional d xy -ba,nds from the zig-zag chains V2-V1'- 
V2-V1' and V2'-Vf-V2'-Vf (Fig. 1), respectively, since 
there are two zig-zag b-chains per unit cell. Each band 
is quarter-filled and one may expect a usual Peierls tran- 
sition to occur due to an appropriate dimerization of the 
structure in b-direction. Although this could explain the 
insulating properties at room temperature, such an expla- 
nation is obviously incorrect, since it would rule out a spin- 
Peierls transition at low temperature. Of course the large 
ratio Udlt X y already excludes the straightforward applica- 
tion of the band picture. Nevertheless for sufficiently large 
t xy one expects that the correlated band picture applies. 
Estimates of the ground state energies for the localized 
regime and the band picture suggest that this is the case 
for t xy > 1.5i a . 

The bands obtained by a real bandstructure calcula- 
tion deviate from this idealistic picture due to small inter- 
chain couplings, which is also reflected in shorter Vl-03- 
V2 bonds (lj]. The full bandstructure calculation predicts 
a metallic state. 

We remark that photoemission experiments in com- 
bination with XAS would provide a sensitive test of the 
present picture, since the lowest unoccupied states are the 
singlet and triplet states, while higher lying 'two-particle' 
states are split by Ud- Optical conductivity measurements 
on the other hand should show a relatively small 'single 
particle' gap of order 2t a - 

In summary we have shown that o/-NaV2 05 is an in- 
sulator due to strong correlations, i.e. large Ud- However 
since the single electron orbitals have bonding character, 
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i.e. involving two V-atoms, the magnetic structure is in- 
duced via 2-electron molecular singlet and triplet states 
which are Heitler-London split. Due to the almost perfect 
cancellation of the triplet- and singlet-interactions in a- 
direction the spin system is essentially one-dimensional. 
We stress that doubly occupied states with energy ~ Ud, 
which usually contribute to the magnetic interaction in 
Mott-Hubbard insulators, have little influence since they 
are at much higher energy. Therefore we suggest the name 
Heitler-London insulator as a more precise characteriza- 
tion for such systems. 

We acknowledge stimulating discussions with J. van 
den Brink, Yu. Grin, A. M. Oles and H.-G. von Schnering. 

References 

1. E. Pytte, Phys. Rev. B 10, 4637 (1974); M. C. Cross and 
D. S. Fisher, Phys. Rev. B 19, 402 (1979). 

2. J. W. Bray et al, Phys. Rev. Lett. 35, 744 (1975); 

3. M. Hase, I. Terasaki and K. Uchinokura, Phys. Rev. Lett. 
70, 3651 (1993). 

4. M. Isobe and Y. Ueda, J. Phys. Soc. Jpn. 65, 1178 (1996). 

5. Y. Fujii, H. Nakao, T. Yoshihama, M. Nishi, K. Nakajima, 
K. Kakurai, M. Isobe. Y. Ueda, and H. Sawa, J. Phys. Soc. 
Jpn. 66, 326 (1997). 

6. M. Weiden, R. Hauptmann, C. Geibel, F. Steglich, M. Fis- 
cher, P. Lemmens, and G. Guntherodt, Z. Phys. B 103, 1 
(1997). 

7. T. Ohama, M. Isobe, H. Yasuoka, and Y. Ueda, J. Phys. 
Soc. Jpn. 66, 545 (1997). 

8. M. Koppen, D. Pankert, R. Hauptmann, M. Lang, M. Wei- 
den, C. Geibel, and F. Steglich, |cond-mat/9710183 . 



P. Horsch and F. Mack: A new view of the electronic structure of the spin-Peierls compound a'-NaV^Os 7 

9. A. Carpy and J. Galy, Acta Crystallgr., Sect B 31, 1481 
(1975). 

10. F. Mila, P. Millet, and J. Bonvoisin, Phys. Rev. B 54, 
11925 (1996). 

11. D. Augier, D. Poilblanc, S. Haas, A. Delia, and E. Dagotto, 
Phys. Rev. B 97, R5732 (1997). 

12. H.-G. von Schnering, Yu. Grin, M. Kaupp, M. Somer, R. 
Kremer, O. Jepsen, T. Chatterji, and M. Weiden, to be pub- 
lished. 

13. M. Fischer, G. Els, P. Lemmens, G. Giintherodt, A. S. 
Mishenko, M. Weiden, R. Hauptmann, C. Geibel, and F. 
Steglich, to be published. 

14. Yu. Grin (private communication). 

15. W. A. Harrison, Electronic Structure and Properties of 
Solids (Freeman, San Francisco, 1980). 

16. H. F. Pen, PhD-thesis Groningen (1997), and private com- 
munication. 

17. Covalency effects lead to diagonal shifts of — 4t a , ~2t a — t^ 
and — 2t„ in the rid — 0, 1 and 2 sectors, respectively, where 
*o ~ tpd/( e d — e P + Ud). Although these shifts are important 
for ground state energies, they do not affect the excitation 
energies considered here. 

18. J. H. Jefferson, H. Eskes, and L .F. Feiner, Phys. Rev. 45, 
7959, (1992). 

19. J. Jaklic and P. Prelovsek, Phys. Rev. B 49, 5065 (1994). 

20. P. W. Anderson, Mat. Res. Bull. 8, 153 (1973), P. Fazekas 
and P. W. Anderson, Phil. Mag. 30, 423 (1974). 

21. S. Eggert, I. Affleck, and M. Takahashi, Phys. Rev. Lett. 
73, 332 (1994). 

22. A. W. Sandvik, R. R. P. Singh, and D. K. Campbell, Phys. 
Rev. B 56, 14510 (1997) . 



